Testing and selection cosmological models with dark energy 
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Abstract. It is described dynamics of a large class of accelerating cosmological models in terms of dynamical systems of 
the Newtonian type. The evolution of the models is reduced to the motion of a particle in a potential well parameterized by 
the scale factor. This potential function can be reconstructed from distant supernovae type la data and many cosmological 
models represented in terms of the potential becomes in a good agreement with current observational data. It is proposed to 
use the information criteria to overcome this degeneracy within a class of A) dark energy models and B) simple models basing 
on modification of the FRW equation. Two class of models can be recommended by the Akaike (AIC) and Schwarz (BIC) 
information criteria: the phantom and ACDM models. 
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oo 1. INTRODUCTION 

m ' 

Recent measurements of distant supernovae (SNIa) [1,2] indicate that our Universe is in an accelerating phase of 
expansion. This has stimulated the search for models explaining this phenomenon. There are two alternative types of 
explanation observed acceleration: 1) the Universe is accelerating due to presence of so called dark energy violating 
the strong energy condition [3], 2) instead of dark energy of unknown form the modification (or generalization) of 
the FRW equation is postulated. While the cosmological constant is still attractive proposition for dark energy, it 
offers only some effective model without understanding of observed value of vacuum energy. This situation open 
the cosmology on different propositions of candidates for dark energy description [4, 5]. All these models can be 
described in terms of a two-dimensional dynamical system of the Newtonian type if it is assumed a general form of the 
equation of state px = w{a)px for dark energy X [6, 7, 8]. From the conservation condition (under condition spatial 
homogeneity and isotropy) 

p = -3H(p+p), (1) 

where p and p are effective energy density and pressure respectively, H = (In a)' is the Hubble function, a is the scale 
factor, a dot denotes the differentiation with respect to the cosmological time. We obtain that p — p(a) or p = p(z), 
where ( 1 + z) = a~ 1 (ao- the present value of the scale factor is chosen as a unity). 

From the Raychaundhuri equation (called also the acceleration equation) we obtain additional equation governed 
evolution of the Universe 

- = - 1 -(p + 3p). (2) 
a o 

One can simply check that equation (2) assumes the form analogous to the form of the Newtonian equation of motion 
a = — if we choose the potential in the form 

pa 2 

V = ~, (3) 
o 

where p plays the role of effective energy density which satisfies the conservation condition. 
Equation (2) admits first integral 

^-+V(a)=E, (4) 

which is called the Friedmann first integral. Finally, due to the particle-like description of the evolution of the 
Universe, which mimics motion of the fictitious particle moving in the potential well. The potential function is 
representing single function of scale factor which determine full evolutional scenarios of the cosmological model. 



Therefore our proposition is using scale factor a and its derivative as independent variables which measure states of 
the system. Moreover in such approach we have additionally Hamiltonian structure and the FRW cosmological system 
is representing in terms of two-dimensional dynamical system. In such formulation space of states is defined in the 
following way: 

e = /(M): ^+V{a) = - k -\, (5) 
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where k is the curvature index, k = 0, ± 1 . £ should be treated as a phase space of the system in which trajectories lies on 

algebraic curves which are determined by the Hamiltonian constrained <ffi = + V(a) = — |. If we reparameterize 
time variable 1 1— > T : \Ho\dt = dx (Hq ^ 0) then £ can be defined equivalently as: 

C = {M):^+VW = -^}, (6) 

where Q^q = / 3 (#o) 2 is the density parameter for curvature fluid, x = ^ is a dimensionless scale factor expressed 
in its present value ao- If we assume that fluid which filled the Universe is composed with noninteracting fluids: dust 
matter, radiation and the cosmological constant then we obtain: 

^f+V(x) = £ f (7) 

where 

Of course all density parameters £2,- o are not independent and satisfy the constraint condition £;=(m. ^ k) Q,0 — 1- 
The original Friedmann first integral (4) can be rewritten to the form 



c/x x 2 
dx 



■ Om,0 (X- 1 - 1 ) + &rfi(x~ 2 - 1 ) + ^A.O^ 2 - 1 ) + 1 • (8) 



"A,0W = 2/- i 7T77 — 71 • W 



In order to survey the set of all its solutions (which have an arbitrary parameter) we pick a pair (i2 m .o,i2 r o) and 
represent the range of each possible solution x(t) corresponding to these data and an initial value £Ia,q, as a horizontal 
segment in a plane (x, i2 \ o)- F° r this purpose it is necessary to mark the boundary of the permitted region (x >Q,x 2 > 
and V(x) < jQ.^ q). The boundary is corresponding to V = j£lk,o an d a — (Qa,o - axis). From V (x) = ^O^ o, 
(x /2 = 0) we obtain: 

x(x - 1 )ft m;0 + (x-l)(x+l )£l rfi - x 2 
x 2 (x+l)(x-l) 

From this boundary curve it is possible to obtain qualitatively a classification of all evolutional scenario in the 
configuration space. Each model corresponding to a maximally extended line segment CI^ q = const which includes 
the present value x = 1. If that segment (£2a,o(*) = cons t level) is bounded on the left by the locus a — 0, the model 
has a big bang. The potential function is given by 

, , l x 2 (x 2 -l)n A ,o+x(l-x)n m; o + (l-x 2 )^ +x 2 
2 x 2 ' 

and all trajectories lies in the classically allowed region V (x) < ji^.o- 

Let us consider now a more general model filled by dark energy satisfying the equation of state px = wx(a)px 
instead of the cosmological constant A. Then we obtain the potential function in the form 

v( x _ 1 x 2 (/(x) - IJOxfl +x(l -x)Cl mfl + (1 -x 2 )n r ,o+x 2 

V(X)- 2 ^ (11) 



where 

w x (aj 



Q. x = £l x ,of(x), f{x) =x 1 exp 



da 



The full classification of dark energy models can be performed in the analogous way like to lambda models from 
consideration of a boundary curve V (x) = 5^,0- 



TABLE 1. Different models explaining acceleration of the Universe in terms of dark energy (A) and modification 
of the FRW equation (B). All models are defined in terms of a single potential function. 



case 



name of model 



form of potential function 



A 1 


models with matter and 


V(x) = 






l-3w^ 






noninteracting fluids 














Pi = WiPi, i=l,...,n 














2 


phantom models 


V{x) = 




l +^x,ox 3 ) 








px = — \px and dust 














3 


models with generalized 
Chaplygin gas 


V{x) = 


1 




A +£l Chfi X 2 


' A 4- l ~ A ' \ ^ 
K A S T t 3(l+c) J 




4 


models with dynamical E.Q.S. 


V(x) = 


4 










PX = {w { )+W X z)Px 















B 5 


RSB models 


V(x) 


= -\ (o. mfl x-i +o. drfi x- 2 + a U) x- 4 + a Afi x 2 ) 


6 


Cardassian models 

P = PmflX + PcfiX 


V(x) 


= -^{a mfl x- 1 +a Cfi x- 3n + 2 ) 


7 


DGP brane models 


V{x) 


= ~h (\J ^mfix- 1 +Q. rcfi X 2 + Q. rc p^ 


8 


bouncing models 


V(x) 


= -\ {a m , x- m + 2 - £2„.ox-"+ 2 + a^x 2 ) 



2. ENSEMBLE OF DARK ENERGY MODELS 



Table 1 shows some examples of different dark energy models appeared in cosmological investigations sources of 
acceleration. In Table 1 we present also potentials function for alternative propositions of explanation acceleration 
without conception of dark energy. For these class of models it is proposed some modification of the Friedmann first 
integral (so called the Cardassian models) or new brane physics is incorporated (Dvali-Gabadadze-Porrati (DGB) 
models or the Randall-Sundrum brane cosmological models). 

Let us introduce the definition of basic notion. By the ensemble of the FRW cosmological models with dark energy 
it is understood a subspace of all dynamical systems on the plane which are represented the FRW dynamical systems of 
a Newtonian type, characterized by the potential function V (a). This space can be naturally equipment in structure of 
the Banach space if we introduce the C 1 metric d(l,2)= max^ e £{|Vi — V2I, \V\ X — ^2*1}, where £ is a closed subspace 
of the configuration space, V\ and V2 are potential functions labeled as 1&2. 

The Bayesian approach to the hypothesis testing was developed by Jeffreys [9, 10] as a program of methodology 
for quantifying the evidence of two competing scientific theories. In our context we test hypotheses about origins of 
the present acceleration of the Universe. In Jeffreys' approach the crucial role plays statistical models which are used 
to represent the probability of the data. The Bayes theorem is used to calculate the posterior probabilities that one of 
the models is correct or to find the most favored model. 

The characteristic feature of Bayesian framework is that it was applied in several fields of sciences. The Bayesian 
methodology offers the possibility of comparison of nonnested models and accounts for uncertainty in the choice of 
models [11]. 



3. IDEA OF MODEL SELECTION 

In the development of cosmology the basic role played an idea of cosmological models together with an idea of 
astronomical tests [12]. The idea of cosmological tests make cosmological models parts of astronomy which offers 
possibility of observationally determining the set of realistic parameters, that can characterize viable models. While we 
can perform estimation of model parameters using a standard minimization procedure based on the likelihood method, 



many different scenarios are still in a good agreement with observational data of SNIa. This problem which appears in 
observational cosmology is called the degeneracy problem. To solve this problem it is required to differentiate between 
different dark energy models. We propose to use the Akaike information criterion (AIC) and Bayesian information 
criterion (BIC). 

For the model selection framework it is required to have data and a set of models and then the model based statistical 
inference. The model selection should be based on a well-justified single (even naive) model or, at least, a simple model 
which suffices for making inferences from the data. In our case the ACDM model plays just the role of such a model 
and application of information criteria seems to be reasonable as a first step toward an extension of the standard 
maximum likelihood method. The model selection should be viewed as a way to obtain model weights, not just a way 
to select only one model (and then ignore that selection occurred). Moreover in the notion of true models do not believe 
information theories because the model by definition is only approximations to unknown physical reality: There is no 
true model of the Universe that perfectly reflect large structure of space-time but some single best model has been 
found. 

In this paper the Bayesian factor is used to select among dark energy models. The Bayesian factor is the posterior 
odds of the one hypothesis when the prior probabilities of two alternative hypotheses are equal. It can be defined as 
the ratio of the marginal likelihoods under two competing models, i.e. the probabilities of the data with all the model 
parameters integrated out. Thus the marginal likelihood is a basic quantity (evidence) needed for this approach [13]. 
We use the Bayesian factor for evidence provided by SNIa data in favor of one cosmological model with dark energy 
(represented by statistical model) as opposed to another. The well known the Bayesian information criterion (BIC) 
gives the approximation to the logarithm of the Bayes factor and can be easily calculated without the requirement of 
prior distribution. It is necessary to apply this method when a very large number of different dark energy models is 
considered and we wish to select only one of them. This problem is solved in the Bayesian framework by choosing the 
model with the highest posterior probability. 

The AIC is defined in the following way [14] 

AIC = -2ln^f + 2d, (12) 

where «£? is the maximum likelihood and d is a number of model parameters. The best model with a parameter set 
providing the preferred fit to the data is that which minimizes the AIC. While there are justification of the AIC in 
information theory and also rigorous statistical foundation for the AIC, it can be also justified as Bayesian using a 
'savvy' prior on models that is a function of a sample size and a number of model parameters. For the AIC we can 
define AAIC, as the difference of the AIC for model ; and the AIC value for the scale model: A, = AIC, — AIQ . The A, 
are easy to interpret and allow a quick 'strength of evidence' comparison and a ranking of candidates for dark energy 
description. The models with A, < 2 have substantial support (evidence), those where 2 < A, < 7 have considerably 
less support, while models having A, > 10 have essentially no support. 

In the Bayesian framework a best model (from the model set {Mi}) is that which has the largest value of probability 
in the light of data (so called posteriori probability): P(M[\D) = ^MM^Mil ; w here P(M;) is a priori probability for 
the model M,, P(D\M{) is the evidence (Ef). We can define the posterior odds for models M, and Mi, which (in the case 
when no model is favored a priori) is reduced to the evidence ratio (so called the Bayes factor - By). 

Schwarz [15] showed that for observations coming from linear exponential family distribution: ln£ , I = InJz? — 
^dlnN + 0(1), where N is number of data points used in the fit (this was extended by Haughton [16] for curved 
exponential family). According to this result Schwarz introduced a criterion for the model selection: the best model is 
that which minimized the BIC, defined as 

BIC = -21n£? + dlnN. (13) 

To compare models M; and Mj one can compute 21nBy = —(BIC, — BIC/) = —ABIC,/ which can be interpret as 
'strength of evidence' against j model: < 21nBy < 2-not worth more than a bare mention, 2 < 21nBy < 6 - positive, 
6 < 21nBy < 10 - strong, and 21nBy > 10 - very strong. 

There are many simulation studies in the statistical literature on either the AIC or BIC alone, or often comparing 
them and making recommendation on which one to use. It should be pointed out that assumptions of using classical 
model selection methods are satisfied, namely: 

1 . there is model-based inference from SNIa data (luminosity distance observable) 

2. there is a set of dark energy models and no certainty about which model should be used in explanation of present 
acceleration 



TABLE 2. Values of the AIC, BIC, AAIC/i and 2\nB u for flat models from 
Table 1. 



Case 








1 In R 


A l(w x = -l) 


179,9 


186,0 


0,0 


0,0 


2 


178,0 


184,1 


-1,9 


-1,9 


3 


181,7 


193,9 


1,8 


7,9 


4 


180,5 


192,7 


0,6 


6,7 


B 5 


182,1 


194,3 


2,2 


8,3 


6 


178,5 


187,7 


-1,4 


1,7 
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1,0 
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FIGURE 1. Value of A) AAIC/j =AIC,-AICi and B) 2 lnB u for models from Table 1. 



3. a data-based choice is made among these competing models (see Table 1) 

Then the inference is made from this one selected model as if it were a priori the only model fit to the data. 

Table 2 gives the value of the AIC and BIC for flat models from Table 1 . It also contains values of AAICn and 
21nZ?i,, where model 1 is our scale model (ACDM). 

From results contained in Table 2 one can conclude that the phantom model (case 2) is the best model from our 
model set, but the ACDM model has still substantial support (evidence). The Akaike information criterion indicates 
that there is only one strong disfavored model (case 8), whereas the Bayesian information criterion denotes that also 
cases 3,4 and 5 are less probable than the ACDM model (see Figure 1). 
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